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Abstract. The Virgo cluster spiral galaxy NGC 4654 is supposed to be a good candidate for ongoing ram pressure 
stripping based on its very asymmetric Hi distribution. However, this galaxy also shows an asymmetric stellar 
distribution. Numerical simulations using ram pressure as the only perturbation can produce a tail structure of 
the gas content, but cannot account for its kinematical structure. It is shown that a strong edge-on stripping 
event can produce an asymmetric stellar distribution up to 800 Myr after the stripping event, i.e. the galaxy’s 
closest passage to the cluster center. Simulations using a gravitational interaction with the companion galaxy 
NGC 4639 can account for the asymmetric stellar distribution of NGC 4654, but cannot reproduce the observed 
extended gas tail. Only a mixed interaction, gravitational and ram pressure, can reproduce all observed properties 
of NGC 4654. It is concluded that NGC 4654 had a tidal interaction ^500 Myr ago and is continuing to experience 
ram pressure. 
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dynamics 


1. Introduction 

Maps of the gas content of spiral galaxies in the Virgo clus¬ 
ter have revealed that the Hi disks of cluster spirals are 
disturbed and considerably reduced (Cayatte et al. 1990, 
1994). Despite their Hi deficiency, these spiral galaxies do 
not show a lack of molecular gas (Kenney & Young 1989). 
These observational results suggest that gas removal due 
to the rapid motion of the galaxy within the intraclus¬ 
ter medium (ICM) (ram-pressure stripping; Gunn & Gott 
1972) is responsible for the Hi deficiency and the disturbed 
gas disks of theses cluster spirals. Nevertheless, it is still 
an open question where and when these galaxies lost their 
atomic gas. 

In order to study in detail how ram pressure acts, we 
are modeling galaxy orbits and use a three-dimensional 
N-body code to simulate the gas kinematics of a spiral 
galaxy falling into the Virgo cluster (Vollmer et al. 2001). 
Ram pressure exerted by the ICM on the ISM of a rapidly 
moving galaxy is explicitly included in this code. 

One candidate for an ISM-ICM interaction in the 
Virgo cluster is NGC 4654. Its physical parameters are 
listed in Tab. |^. Phookun & Mundy (1995) observed 
this galaxy in the 21-cm-line with a high sensitivity 
(^ 10^® cm“^) at the VLA. They found a very asym¬ 
metric Hi distribution with a compressed edge on one side 
and a long tenuous tail on the other side. They argued 
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Table 1. Physical Parameters of NGC 4654 


Other names 


a ( 2000 )’^ 

5 (2000)’" 

Morphological type*" 

Distance to the cluster center (°) 
Optical diameter Djs*" (^ 

B O a 
T 

Systemic heliocentric velocity*" (hms”*^) 
Distance D (Mpc) 

Vrotmax (kms“^) 

PA 

Inclination angle 
HI dehciency*^ 


UGC 7902 
VCG 1987 
GGCG 071 - 019 
12‘’43“56.6= 
13°07'36" 
SBcd 
3.4 
4.9 
10.75 
1035±3 
17 

lOO*", 139” 
121° 128° ° 
49° 51° ° 

0.17±0.2 


*" RC3, de Vaucouleurs et al. (1991) 
Guharthakurta et al. (1988) 

° Sperandio et al. (1995) 

Cayatte et al. (1994) 


that this asymmetry is due to ram pressure stripping that 
pushes the ISM beyond the galaxy’s optical radius. 

However, NGC 4654 does not only show an asymmetric 
gas disk, but also has an asymmetric stellar distribution. 
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This might be due to a gravitational interaction with the 
nearby spiral galaxy NGC 4639. 

In order to investigate the possibilities of tidal or 
ISM-ISM interaction, the code presented in Vollmer et 
al. (2001) has been extended and contains now a non- 
collisional component. This gives the opportunity to make 
simulations of a gravitational interaction with or without 
additional ram pressure. 

In this article I compare snapshots of four different 
simulations with Hi observations of NGC 4654: 

1. ram pressure, analytically given gravitational potential 
of the galaxy; 

2. ram pressure, galactic halo, bulge, and disk simulated 
by non-collisional particles; 

3. gravitational interaction; 

4. gravitational interaction and ram pressure. 

The comparison of the gas distribution and velocity held 
permits to discriminate between the different scenarios. 

The plan of this article is the following: the numerical 
code is described in Sect. 2. The simulations using ram 
pressure as the only perturbation and their comparison 
with observations are presented in Sect. 4. In Sect. 5 a 
scenario of active stripping is discussed. The simulations 
using a gravitational interaction with and without ram 
pressure and their comparison with observations are pre¬ 
sented in Sect. 6, followed by a discussion of the results of 
all simulations in Sect. 7. A summary and conclusions are 
given in Sect. 8. 

2. The model 

2.1. The collisional component 

Since the model is described in detail in Vollmer et al. 
(2001), I will only briefly summarize its main features. The 
particles represent gas cloud complexes which are evolv¬ 
ing in an analytically given gravitational potential of the 
galaxy. 

10 000 particles of different masses are rotating within 
this gravitational potential. The total gas mass is Mg°l = 
5.810® Mq, which is ~30% larger than the observed gas 
mass (Huchtmeier & Richter 1989). To each particle a ra¬ 
dius is attributed depending on its mass. During the disk 
evolution the particles can have inelastic collisions, the 
outcome of which (coalescence, mass exchange, or frag¬ 
mentation) is simplified following Wiegel (1994). This re¬ 
sults in an effective gas viscosity in the disk. 

As the galaxy moves through the ICM, its clouds are 
accelerated by ram pressure. Within the galaxy’s iner¬ 
tial system its clouds are exposed to a wind coming from 
the opposite direction of the galaxy’s motion through the 
IGM. The temporal ram pressure profile has the form of a 
Lorentzian, which is realistic for galaxies on highly eccen¬ 
tric orbits within the Virgo cluster (Vollmer et al. 2001). 
The effect of ram pressure on the clouds is simulated by 
an additional force on the clouds in the wind direction. 
Only clouds which are not protected against the wind by 
other clouds are affected. 


Table 2. Number of particles N, particle mass M, and 
smoothing length I for the different galactic components. 


component 

N 

M (Mo) 

1 (pc) 

halo 

16384 

9.2 10“ 

1200 

bulge 

16384 

3.2 10® 

180 

disk 

16384 

1.610® 

240 

companion 

11000 

9.2 10® 

1200 


2.2. The analytical gravitational potential 

The fixed, analytical gravitational potential consists of 
three parts: the dark matter halo, the stellar bulge, and 
the disk potential (Allen & Santillan 1991). Using their 
dehnitions, the model parameters are: (i) halo: 03=12 kpc 
(effective radius), M 3 = 8.610^® M© (halo mass), (ii) 
bulge: 61=390 pc (bulge core radius). Mi = 5.610® Mq 
( bulge mass), (iii) disk: 02=2.7 kpc (disk scale length), 
62=250 pc (disk scale height), M 2 = 2.610^® M© (mass 
of the stellar disk). The resulting velocity held has a con¬ 
stant rotation curve of Urot "^150 kms“^, which is consis¬ 
tent with the rotation curves derived by Guharthakurta 
et al. (1988) and Sperandio et al. (1995). 

2.3. The non-collisional component 

The non-collisional component consists of 49125 parti¬ 
cles, which simulate the galactic halo, bulge, and disk. 
The characteristics of the different galactic components 
are shown in Tab. ^ The particle trajectories are inte¬ 
grated using an adaptive timestep for each particle. This 
method is described in Springel et al. (2001). The follow¬ 
ing criterion for an individual timestep is applied: 


where ai is the acceleration of the particle i. The minimum 
of all ti dehnes the global timestep used for the Burlisch- 
Stoer integrator that integrates the collisional component. 

The setup of the initial conditions was made by the 
program described in Hernquist (1993) with the following 
parameters: 

— halo: mass M = 1.510^^ Mq, core radius rc=3 kpc, 
tidal radius rt=30 kpc, cutoff radius i?=90 kpc, 

— bulge: mass M = 5.2 10® Mq, scale length Z=0.3 kpc, 
cutoff radius i?=30 kpc, 

— disk: mass M = 2.610^® Mq, scale length Z=3 kpc, 
cutoff radius i?=45 kpc, disk thickness 2;o=600 pc. 

10 000 collisional particles with a 1/R column density 
profile were added and the system was evolved during 
2 Gyr in order to obtain a relaxed system. At the end of 
the simulation the difference between the total energy and 
the total angular momentum of the system was smaller 
than 0.5% of their initial value. During the last Gyr of 
this simulation the disk properties (disk height, surface 






B. Vollmer: NGC 4654 


3 


density profile, density profile) did not change within 5%. 
The final state of this simulation was used as the initial 
state for the here presented simulations. 

The companion galaxy with a mass of 10^^ Mq is sim¬ 
ulated by 11000 particles forming a Plummer sphere with 
with a core radius of 2 kpc. 

3. The simulations 

I have made four different simulations using different inter¬ 
actions and different methods to model the galactic com¬ 
ponents (dark halo, stars, and gas). 

1. Time dependent ram pressure stripping, no gravita¬ 
tional interaction. The gravitational potential of the 
galaxy is fixed and analytically given. 

2. Time dependent ram pressure stripping, no gravita¬ 
tional interaction. The gravitational potential of the 
galaxy is simulated by non-collisional particles. 

3. Gravitational interaction, no ram pressure stripping. 
The gravitational potential of the galaxy is simulated 
by non-collisional particles. 

4. Gravitational interaction and constant ram pressure. 
The gravitational potential of the galaxy is simulated 
by non-collisional particles. 

4. Time dependent ram pressure stripping 

In Vollmer et al. (2001) we simulated different galaxy or¬ 
bits with different inclination angles i between the or¬ 
bital and the disk plane. Phookun & Mundy (1995) sug¬ 
gested that NGG 4654 is moving edge-on through the 
IGM. Indeed, only simulations with i < 20° show an ex¬ 
tended tail several 10® yr after the galaxy’s closest pas¬ 
sage to the cluster center. Furthermore, the simulation 
has to reproduce the observed Hi deficiency of NGC 4654. 
With i = 0° (edge-on stripping), a maximum ram pres¬ 
sure > 5000 cm“®(kms“^)^ leads to a final Hi de¬ 
ficiency DEF > 0.05. It turned out that the best corre¬ 
spondence between the observed and the model Hi surface 
density exists for Pram = 5000 cm“®(kms“^)^. In Vollmer 
et al. (2001) we have shown that the highest realistic max¬ 
imum ram pressure is 10000 cm“®(kms“^)^. In 

order to obtain higher values the galaxy would have to ap¬ 
proach M87 to distances where a damaging gravitational 
interaction is ineluctable. Thus, values of the maximum 
ram pressure in the range 5000 cm“®(kms“^)^ < Pram < 
10000 cm“®(kms“^)^ and the inclination angles in the 
range of 0° < i < 10° can reproduce the Hi observations. 
Since the observed Hi deficiency has an error of ±0.2, the 
final model deficiency of DEF=Q.05 is consistent with the 
observed value {DEF=0.17), i.e. no Hi deficiency. The 
best fit model parameters are given in Tab. ^ 

4.1. The projection on the sky 

In order to project the model gas distribution on the 
sky, three angles are needed: the position angle, the in¬ 
clination angle, and the azimuthal angle within the plane 


Table 3. Model parameters for time dependent ram pres¬ 
sure stripping 


maximum rampressure (cm“®(kms~^)^) 5000 

inclination angle between orbital and disk plane 0° 
final HI deficiency 0.1 


of the galactic disk. The position angle can be deter¬ 
mined observationally. If one assumes that the spiral arms 
of NGC 4654 are trailing, its velocity field (Phookun & 
Mundy 1995) places its northern edge in front of the 
galaxy center, i.e. the sign of the inclination angle i is 
known. 

In general, galaxies have trailing spiral arms. Galaxies 
with leading spiral arms are assumed to have undergone a 
plunging retrograde encounter with another galaxy (see, 
e.g., Thomasson et al. 1989 or Byrd et al. 1993), which 
should not be the case for NGC 4654. Since NGC 4654 
harbors a bar and two normal spiral arms, there is no 
apparent reason why its spiral arms should be leading. It 
will be shown in Sect. || that the assumption of trailing 
spirals is important for the model. 

With a given position angle (PA) and a given incli¬ 
nation angle *, only the azimuthal (a) in the galactic 
plane can be varied. I assume that the galaxy orbit within 
the cluster is approximately linear in space (see Vollmer 
et al. 2001). If NGC 4654 is emerging from the cluster 
core, it must be located in front of M87, because its ra¬ 
dial velocity with respect to the cluster mean is negative 
{Av ~ —100 kms“^). Furthermore, the component of its 
three-dimensional velocity vector parallel to the right as¬ 
cension axis must be negative (the galaxy is moving to the 
east). 

In Fig. I^the three components of the galaxy’s velocity 
vector are shown as a function of a. is the com¬ 
ponent parallel to the axis of right ascension (positive is 
to negative right ascensions), is that parallel to the 
axis of declination, and is that parallel to the line of 
sight. The constraints < 0 kms“^ and < 0 kms“^ 
imply 130° < a < 220°. It turned out that a = 180° leads 
to the best fit to observations. This angle is marked by an 
arrow in Fig. |^. 

4.2. The simulations 

Fig. II and ^ show the evolution of the galaxy’s ISM after 
its closest passage to the cluster center for the model with 
a fixed gravitational potential and a galaxy model includ¬ 
ing a non-collisional component respectively. The position 
and inclination angle of NGG 4654 are used. The time t 
can be found above each frame. The closest passage to 
the cluster center corresponds to t=0 yr. The arrow indi¬ 
cates the direction of the wind, i.e. it is opposite to the 
galaxy’s motion within the cluster. The length of the ar¬ 
rows is proportional to the ram pressure Pram = PiCM'^gao 
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Fig. 1. The three components of the galaxy’s velocity vec¬ 
tor as a function of the azimuthal angle a. Solid: 
dotted: dashed: The chosen angle is marked by 

an arrow, is the component parallel to the axis of right 
ascension (positive is to negative right ascensions), is 
that parallel to the axis of declination, and is that 
parallel to the line of sight. 

where picM is the ICM density and Ugai is the galaxy ve¬ 
locity with respect to the cluster mean. The northern edge 
of the galaxy is in front of the galaxy center. The galaxy 
is rotating counter-clockwise. Since both simulations are 
very similar, I will describe in the following the main fea¬ 
tures for both simulations: at approximately maximum 
ram pressure (t = 0) an overdensity builds up in the di¬ 
rection of the galaxy’s motion, where the gas has been 
pushed to smaller radii. At t ~ 100 Myr the wind has 
driven out the gas to the north where the vector of the 
galactic rotation velocity is parallel to the wind direction. 
In Vollmer et al. (2001) this is called the accelerated arm. 
A weak decelerated arm can also be seen in the south. 
At t ^ 200 Myr ram pressure has already ceased com¬ 
pletely and the evolution of the galaxy is entirely due 
to rotation and re-accretion of the gas, which has not 
been accelerated to the escape velocity (uesc 
Due to rotation, the accelerated arm moves to the north 
and a part of its material falls back onto the galaxy. At 
t > 350 Myr an asymmetric shell structure is building up 
in the south, which is most prominent in the south-east. 
This south-eastern part expands and forms an extended 
tail at t ^ 700 Myr after the galaxy has passed the cluster 
center. 

The main difference of the gas dynamics between the 
two models is the that the gravitational potential of model 
galaxy including a non-collisional component can change 
with time. The ISM-ICM interaction pushes the gas to 
smaller galactic radii and heats it. This process also heats 
the stellar disk. At maximum ram pressure the surface 
density in the north-east increase by a factor ^2 and 
~2.5 10® Mq of gas, which represents ~10% of the total 
enclosed mass within 5 kpc, are displaced to the south¬ 



-40 -20 0 20 40 -40 -20 0 20 40 

R* (kpf] R4 (kpc) 


Fig. 2. Snapshots of the evolution of the galaxy’s ISM 
with a fixed gravitational potential. The cc-axis corre¬ 
sponds to the right ascension. The y-axis to declination. 
The elapsed time is indicated at the top of each panel. 
The position and inclination angle of NGC 4654 are used. 
The galaxy rotates counter-clockwise. It is moving to the 
north-east, i.e. the wind is coming from the north-east 
(indicated by the arrows). The length of the arrow is pro¬ 
portional to the ram pressure (picM'*^gai)- 

west. In addition, this gravitational perturbation triggers 
the formation of an asymmetry of the stellar distribution, 
which can also heat the stellar disk. Thus, the length scale 
of its gravitational potential in z direction increases, i.e. 
the motion of the gas is less two dimensional than in the 
case of the fixed gravitational potential. The result is that 
the asymmetries of the gas distribution are smoothed out 
in the case of the model including a non-collisional com¬ 
ponent. In a real galaxy the gas might cool through shocks 
during the interaction. Moreover, the disk heating through 
the the heated gas might be a numerical artifact due to 
the relatively small number of particles used in the simula¬ 
tion. It is not clear if the gravitational heating is sufficient 
to produce the large scale height observed in the model. 

There is another big difference between the two mod¬ 
els: the disk and bulge stars react to the change of the 
gravitational potential when the gas is pushed to small 
galactic radii. This produces a gravitational shock, i.e. the 
stellar surface density also shows an asymmetry with a 
delay of ^400 Myr. This asymmetric distribution rotates 
slowly and smears out within a few rotation times. Fig. ^ 
shows the distribution of disk stars seen face-on and pro- 
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Fig. 3. Snapshots of the evolution of the galaxy’s ISM 
for a galaxy model including a non-collisional compo¬ 
nent. The x-axis corresponds to the right ascension, The 
j/-axis to declination. The elapsed time is indicated at the 
top of each panel. The position and inclination angle of 
NGC 4654 are used. The galaxy rotates counter-clockwise. 
It is moving to the north-east, i.e. the wind is coming from 
the north-east (indicated by the arrows). The length of the 
arrow is proportional to the ram pressure (picM'^gai)- 


jected on the sky with the parameters of NGC 4654 at 
t=800 Myr, i.e. the snapshot that I want to compare with 
observations. The bottom panel of Fig. ^ shows the K 
band image of Boselli et al. (1997). Seen face-on, the stel¬ 
lar surface density of the outer disk (i? > 6 kpc) is higher 
in the north than in the south. The stellar density in the 
inner disk is high with a strong decline at i? ^ 6 kpc. This 
translates into an elliptical projected stellar distribution 
with the galaxy center in the south-eastern focal point. 
This resembles the overall observed distribution of the old 
stellar component in K band (Fig. || bottom panel). Thus 
the observed overall asymmetry of the stellar distribution 
can be reproduced by a gravitational shock due to edge- 
on ram pressure stripping, but it is not possible to create 
a bar nor spiral arms in this way. 

Fig. I shows the stellar distribution and its velocity 
field at t=800 Myr projected on the sky with the parame¬ 
ters of NGC 4654. The velocity field shows clearly a steep¬ 
ening in the north-west, whereas the southern part of the 
galaxy has the velocity field of a rotation curve that flat¬ 
tens in the outer disk. 





Fig. 4. Top panel: model distribution of the disk stars at 
t=800 Myr seen face-on. Middle panel: model distribution 
of the disk stars at t=800 Myr projected on the sky with 
the PA and i of NGC 4654. Bottom panel: K band image 
of NGC 4654 (Boselli et al. 1997). 




Fig. 5. Left panel: model stellar distribution at 
t=800 Myr. Right panel: model stellar velocity field. Both 
distributions are projected on the sky using PA and i of 
NGC 4654. 


4.3. Comparison between observations and simulations 


In this Section the gas distribution and the velocity held of 
NGC 4654 are directly compared to the last model snap¬ 
shots of Fig. ^ and i.e. at t=800 Myr. 

Fig. ^ shows the Hi distribution of NGC 4654 (bottom 
panel) together with the model gas distributions convolved 
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Fig. 6. Top panel: model distribution of the Hi gas using 
a fixed gravitational potential. The contours correspond 
to (1, 5 ,9 ,13, 17, 21, 25)xl0^® cm“^. Middle panel: 
distribution of the Hi gas using the model including a 
non-collisional component with the same contour levels. 
Bottom panel: Hi distribution of NGC 4654 (Phookun & 
Mundy 1995) with the same contour levels. 



to the beamsize of the 21-cm-line observations (25"). The 
top panel shows the gas distribution of the model with 
a fixed gravitational potential, the middle panel shows 
that of the model including a non-collisional component. 
The main characteristics of the observed gas distribution 


(Fig. ^ bottom panel) are (i) the prominent maximum at 
the north-western edge of the disk, (ii) the overall asym¬ 
metry within the disk with more gas in the south-east, 
(iii) the linear extended tail to the south-east. The model 
using a fixed gravitational potential (Fig. ^ top panel) 
can reproduce feature (i) and (ii). It also shows an ex¬ 
tended tail to the south-east, but it has a higher column 
density, is curved, and is much less extended than the 
observed tail. The model including a non-collisional com¬ 
ponent (Fig. ^ middle panel) fails in reproducing feature 
(i) and (ii), but can reproduce feature (iii), i.e. the ex¬ 
tended, low column density, south-eastern tail. However, 
the model tail is slightly curved. The low column density 
of this tail is due to the tidal heati ng o f the disk during 
the ICM-ISM interaction (see Sect. 4.2). 


The Hi velocity fields of the two models and the ob¬ 
served velocity field can be seen in Fig. Q. The main char¬ 
acteristics of the observed Hi velocity field (Fig. ^ bottom 
panel) are (i) the plateau of constant rotation velocity in 
the south-east, (ii) the steepening of the rotation curve 
in the north west, (iii) the structure of the velocity in 
the extended south-eastern tail. The model using a fixed 
gravitational potential (Fig. top panel) shows the fea¬ 
tures (i) and (ii), but fails in reproducing feature (iii). The 
model including a non-collisional component (Fig. mid¬ 
dle panel) shows a less clear plateau, but feature (ii) is also 
visible. It also fails in reproducing the observed velocity 
field in the extended south-eastern tail. 


In order to compare the characteristic structure of the 
Hi velocity fields, the position-velocity diagrams along (i) 
the major axi^ (ii) the tail, and (iii) the minor axis are 
shown in Fig. |8| g, and respectively. The model cut for 
a fixed gravitational potential along the major axis (Fig. |[ 
top panel) nicely shows the observed plateau of constant 
rotation velocity. At the opposite side, the main compo¬ 
nent is spatially reduced with respect to the side with 
positive offsets and shows a constant rotation. Thus, this 
simulation is able to reproduce the south-eastern part of 
the position-velocity diagram, but not the north-western 
part. 


The cut for the model including a non-collisional com¬ 
ponent (Fig. H, middle panel) also shows a plateau at the 
south-eastern side (negative offsets) that is spatially re¬ 
duced having only one third of the extent of the observed 
plateau. The linear rise at the north-western side is re¬ 
produced. At the north-western edge the rotation curve 
is falling again producing closed contours in the veloc¬ 
ity field. This is due to the velocity dispersion of the gas 
within the disk, which is higher than that of the observed 
galaxy, because there is no explicit cooling mechanism in¬ 
cluded in the model. 

The model cuts through the tail (Fig. top and middle 
panel) are very similar for both models. As already seen 
above, these model cuts differ from the observed one in the 
sense that the velocity structure of the model tail is almost 
linear, whereas that of the observed tail is significantly 
curved. 
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Fig. 7. Top panel: model Hi velocity field using a fixed 
gravitational potential. The contours, in kms“^, are from 
900 (north-west) to 1160 (south-east) in steps of 20 for 
both velocity fields. Middle panel: gas velocity field us¬ 
ing the model including a non-collisional component with 
the same contour levels. Bottom panel : Hi velocity field 
of NGC 4654 (Phookun & Mundy 1995) with the same 
contour levels. 


The model cuts along the minor axis (Fig. shows 
qualitatively the same deviations from circular rotation as 
the cut through the observed data cube. 




Offset (orcseconds) 



Fig. 8 . Position-velocity cut through the Hi cube along 
the major axis. Top panel: model using a fixed gravita¬ 
tional potential. Middle panel: model including a non- 
collisional component. Bottom: Hi (Phookun & Mundy 
1995). 


5. The active stripping scenario 


The simulations of Sect. 4.2 show that ram pressure pushes 
gas beyond the optical radius of the galaxy at t > 0 Myr. A 
tail-like structure is built up at even later times, clearly 
after the galaxy’s closest passage to the cluster center. 
This result rules out the scenario of active ram pressure 
stripping. 

Nevertheless, if one assumes that the efficiency of ram 
pressure to push the atomic gas has been greatly un¬ 
derestimated and that the galaxy is moving through a 
region of enhanced ICM density, it is possible to com¬ 
pare another simulation snapshot with the observations, 
where ram pressure stripping is still active. I have cho¬ 
sen the following simulation parameters to best repro¬ 
duce the observed Hi distribution with the extended tail: 
maximum ram pressure Pram=2000 cm“^(kms“^)^, and 
i = Q° (edge-on stripping). The timestep of the snapshot 
is t ^ 100 Myr, i.e. ram pressure is still active. The gas dis¬ 
tribution and velocity field of this snapshot can be seen in 
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Fig. 9. Position-velocity cut through the Hi cube along 
the tail. Top panel: model using a fixed gravitational po¬ 
tential. Middle panel: model including a non-collisional 
component. Bottom: Hi (Phookun & Mundy 1995). 


Fig. 10. Position-velocity cut through the Hi cube along 
the minor axis. Top panel: model using a fixed gravita¬ 
tional potential. Middle panel: model including a non- 
collisional component. Bottom: Hi (Phookun & Mundy 
1995). 


Fig. H. The tail represents the accelerated arm (cf. Fig. g). 
The assumption that the spiral arms are trailing places 
it imperatively in the north (the galaxy rotates counter¬ 
clockwise) . The column density is uniform over the whole 
length of the tail, which is not observed. Concerning the 
velocity field, the southern plateau of constant velocity is 
inclined with respect to the optical major axis. Thus, this 
snapshot clearly does not fit the Hi observations appro¬ 
priately. 


If one goes still further and drops the assumption that 
the spiral arms of NGC 4654 are trailing (Sect. ^ ), the 
accelerated arm in Fig. would be located in the south 
as observed. In this case, two discrepancies between the 
model and the observations persist: (i) the uniform col¬ 
umn density of the model tail and (ii) the inclination of 
the southern plateau of constant rotation velocity of the 
model. 


6. Gravitational interaction 

6.1. Gravitational interaction alone 

In this simulation we assume that NGC 4654 has been 
gravitationally perturbed by its companion NGC 4639. 
Since we need a close encounter in order to generate the 
observed asymmetry of the stellar disk of NGC 4654 with¬ 
out disrupting it, I choose a retrograde encounter. The 
perturbing galaxy is modeled as a plummer sphere with a 
mass of 10^^ Mq and a length scale of 2 kpc. The impact 
parameter is 21 kpc and the maximum relative velocity 
is 350 kms“^. The minimum distance in the plane of the 
disk of NGC 4654 is 18 kpc. I assume here that the incli¬ 
nation angle of the disk of NGC 4639 is such that there is 
no ISM-ISM interaction. 

I define the following coordinate system: x — y plane: 
disk plane, z axis: perpendicular to the disk plane. The 
origin is the center of NGC4654. In this system the ini¬ 
tial conditions for NGC4639 are: a::=-40 kpc, y=-16 kpc, 
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Fig. 11. Top: gas distribution of a simulation snapshot 
at t ~ 100 Myr with a maximum ram pressure of 
2000 cm“^(kms“^)^ and i = 0° (edge-on stripping). The 
contours correspond to (1, 5 ,9 ,13j 17, 21, 25)xl0^® cm“^. 
Bottom: corresponding velocity field. The contours, in 
kms“^, are from 940 (north-west) to 1140 (south-east) 
in steps of 20. 


z=70 kpc, Vx=0 kms“^, Vy=130 kms“^, u^=-130 kms“^. 
It is marked as a cross on the snapshots. 

At t=l Gyr, the projected distance and radial velocity 
difference with respect to NGC 4654 using the PA and 
i of NGC 4654 are: i=60 kpc, y=55 kpc, Ur=70 kms“^. 
This compares to the observed projected position and ra¬ 
dial velocity with respect to NGC 4654 of iobs=80 kpc, 
2/obs=39 kpc, 'y°^®=-35 kms“^. Both, positions and radial 
velocity, are close enough to suppose that the simulation 
should show the main characteristics of a possible gravi¬ 
tational interaction between NGC 4654 and NGC 4639. 

Fig. 1^ shows eight timesteps of this simulation. The 
timestep of closest encounter is t=460 Myr. After the pas¬ 
sage of NGC 4639, NGC 4654 begins to form a large 
bar with two asymmetric spiral arm, the south-eastern 
being more extended than the north-western arm. After 
^350 Myr these large scale spirals begin to wind up. A 


u 
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Fig. 13. Close-up of the stellar disk of the final snapshot 
after 1 Gyr (see Fig 


small, negligible number of stellar particles is initially lo¬ 
cated outside the gas disk. They are responsible for the 
low density outer spiral arms for t > 700 Myr that are not 
present in the gas distribution. On the other hand, the gas 
distribution has a 1 /R profile, whereas the stellar distribu¬ 
tion has an exponential profile. This causes the difference 
between the high column density structure of the gaseous 
and stellar distribution. In the inner part of the galaxy, i.e. 
within its optical radius, a bar with two asymmetric spiral 
arms survive until t=l Gyr, i.e. ~550 Myr after the closest 
encounter. The north-western spiral arm is much tighter 
than the south-western one as observed for NGC 4654 in 
the optical. Fig. shows a closer view of the stellar disk 
at t=l Gyr. The bar and the asymmetric spiral arms can 
be recognized. This morphology resembles closely that ob¬ 
served in the K band (Fig. jj, bottom panel). 

6.2. Gravitational interaction and ram pressure 

As a next step, ram pressure is included in the model: 
Pram = Wp’T-iCM'Cgai, where TOp is the proton mass, tt-icm 
the density of the intracluster medium, and Ugai the veloc¬ 
ity of the galaxy within the cluster. I have chosen it to be 
constant in order to keep the model as simple as possible. 
The values for the density of the intracluster medium and 
the galaxy’s velocity are typical for a galaxy located at the 
projected distance of NGC 4654, i.e. njcM = 2 10“^ cm“^ 
and z;gai=1000 kms“^. The direction of the galaxy’s mo¬ 
tion is parallel to its major axis and is pointing to the 
north-west. The perturbing mass simulating NGC 4639 
has the same orbits as described in Sect. 

Fig. 11 shows eight timesteps of this simulation. 
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Fig. 12. Snapshots of the tidal interaction between NGC 4654 and NGC 4639. The left panels show the stellar disk, 
the panels show the gas. The timesteps are marked on the upper right of each panel showing the gas. NGC 4639 is 
marked as a cross. 
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Fig. 14. Snapshots of the tidal interaction between NGC 4654 and NGC 4639 including a constant ram pressure. The 
left panels show the stellar disk, the panels show the gas. The timesteps are marked on the upper right of each panel 
showing the gas. NGC 4639 is marked as a cross. 
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Already at t=300 Myr ram pressure pushes the gas 
to the south-east of the galaxy that has been forced to 
larger galactic radii due to the tidal interaction with the 
perturbing mass. It forms a gaseous spiral arm at the west¬ 
ern side of the galaxy at t=400 Myr which then moves to 
the north. It arrives at the northern edge of the galaxy at 
t=700 Myr and proceeds further to the south-east, where 
it forms an extended tail structure at t=l Gyr, which 
is corresponds to the observed tail of NGC 4654. As ex¬ 
pected, the stellar distributions at this final timestep for 
the simulation with and without ram pressure (Fig. ^ and 
Fig. 1^ are the same since the gas displacement is not im¬ 
portant enough t o tr igger a strong stellar asymmetry as 
observed in Sect. 4.3. 
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6.3. Comparison between observations and simulations 


In this Section the gas distribution and the velocity field of 
NGC 4654 are directly compared to the last model snap¬ 
shots of Fig. |l^ and |^, i.e. at t=l Gyr. 

Fig. 15 shows the Hi distribution of NGC 4654 (bottom 
panel) together with the model gas distributions convolved 
to the beamsize of the 21-cm-line observations (25"). The 
top panel shows the gas distribution of the model of a grav¬ 
itational interaction alone, the middle panel shows that a 
gravitational interaction together with constant ram pres¬ 
sure. For both models the asymmetry of the gas distribu¬ 
tion within the optical radius along the major axis can be 
reproduced. The observed local maximum in the north¬ 
west as well as the extended region in the south-east are 
clearly visible in both model snapshots. The simulations 
using only ram pressure as perturbation could not repro¬ 
duce these features (see Fig. ||). 

The very extended south-eastern tail can only be re¬ 
produced by the model including ram pressure. 

The velocity fields of the two models and the observed 
velocity field can be seen in Fig. E 

In order to compare the characteristic structure of the 
velocity fields, the position-velocity diagrams along (i) the 
major axis, (ii) the tail, and (iii) the minor axis are shown 
in Fig. 1^, |I^ and respectively. The south-eastern side 
(negative offsets) of the model cuts along the major axis 
are very similar for both models, i.e. the dynamics of this 
part of the disk are not affected by ram pressure. The 
model cuts show a maximum at an offset of ~-50”. The 
rotation curve then decreases slightly. The observed ro¬ 
tation curve does not show a maximum. It reaches the 
plateau of constant rotation velocity at an offset of -40”. 
The north-western side (positive offsets) of the galaxy 
shows a slightly rising rotation curve for the model includ¬ 
ing only a gravitational interaction. In contrast, the model 
of the mixed interaction (gravitational and ram pressure) 
leads to a linear rising rotation curve in the north-west. 
Again, at the north-western edge the rotation curve is 
falling due to the velocity dispersion of the gas within the 
disk, which is higher than that of the observed galaxy, be¬ 
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Fig. 15. Top panel: model distribution of the Hi gas for 
a model of a gravitational interaction alone. The contours 
correspond to (1, 5 ,9 ,13, 17, 21, 25)xl0^® cm“^. Middle 
panel: distribution of the Hi gas of the model of a gravi¬ 
tational interaction together with constant ram pressure. 
Bottom panel: Hi distribution of NGC 4654 (Phookun & 
Mundy 1995) with the same contour levels. 



cause there is no explicit cooling mechanism included in 
the model (cf. Sect. |^ ). 

Since the model of a mixed interaction and observa¬ 
tions show a constant rotation curve at one side and a 
rising rotation curve at the other side, I conclude that the 
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main characteristics of the observed rotation curve are re¬ 
produced by this model. 

The p-V cuts along the tail reflect the structure of 
the gas distributions (Fig. |^). The gas distribution of 
the model without ram pressure extends further to the 
north-west (positive offsets) than the one including ram 
pressure. At the opposite side the situation is reversed. 
The model including ram pressure reproduces nicely the 
observed p-V cut along the tail. The p-V cuts of both 
models along the minor axis are very similar. They show 
a smaller but visible degree of asymmetry compared to 
the observed asymmetry. 

7. Discussion 

The main differences between the two models using ram 
pressure as the only perturbation are that (i) the disk is 
heated during the ICM-ISM interaction and (ii) the stel¬ 
lar distribution reacts to the heavily perturbed gas dis¬ 
tribution. Whereas (ii) is a robust result, the disk heat¬ 
ing might be artificial, because of the discreteness of the 
model and because there is no gas cooling included in the 
model. Thus, the decrease of the gas column density in 
the extended tail has to be interpreted with caution. If 
one wants to compare this kind of models to observations 
of other galaxies, both models have to be taken into ac¬ 
count. The situation is still more complicated by the fact 
that at the outer regions of the gas distribution evapo¬ 
ration effects might play an important role that can in 
principle lower the gas column density. 

The comparison between the models and observations 
shows important differences. Whereas the model using 
ram pressure as the only perturbation can explain the ob¬ 
served asymmetric rotation curve and the extended, low 
surface density tail of NGC 4654, it cannot reproduce the 
observed velocity structure of the extended, south-eastern 
tail. Nevertheless, it is worth noting that a past, strong 
ram pressure event can produce a tail-like structure even 
800 Myr after the galaxy’s closest passage to the cluster 
center. Since we observe a reaction of the stellar disk to 
the strongly perturbed gas distribution, the perturbation 
might be longer-lived than expected if one takes only the 
timescale of gas diffusion into account. Moreover, it is not 
possible to form a bar and asymmetric spirals as observed 
with ram pressure as the only perturbation. 

The model using a gravitational interaction is able to 
reproduce the observed stellar distribution of NGC 4654. 
The chosen orbit of the perturbing galaxy, NGC 4639, is 
not unique. It has been chosen such that (i) the position of 
the perturber and its velocity are close to those observed 
for NGC 4639 and (ii) to reproduce the asymmetry of the 
stellar distribution. The aim of this work is not to explain 
all details of the observations, but to reproduce the main 
characteristics of the observations. It becomes clear that 
only the combination of a gravitational interaction and 
ram pressure can account for all characteristics. I want to 
stress here that it is important to compare the gas distri¬ 
bution and the gas kinematics. In the case of NGC 4645 it 


is the gas kinematics that discriminate between the mod¬ 
els. Since the tidal interaction pulls gas to larger galactic 
radii, only a small amount of ram pressure is needed to 
produce the very extended gas tail. 

Despite the fact that only the mixed perturbation 
(gravitational and ram pressure) can explain all character¬ 
istics of the observations, there are two potential problems 
with this model: 

— is an impact parameter between the two galaxies of 
~20 kpc realistic? In this case the disk of NGC 4639 
has to be almost parallel to that of NGC 4654 during 
the encounter to avoid an ISM-ISM interaction, which 
would complicate the scenario. This can be verified by 
making a model using two galaxy models that include 
a collisional and a non-collisional component. This in¬ 
troduces a further open parameter, i.e. the inclination 
angle of disk of NGC 4639 with respect to the disk of 
NGC 4654. 

— In the restframe of NGC 4654, NGC 4639 has a ve¬ 
locity in the plane of the sky of fN 4639 ~ 145 kms“^ 
towards the north-west. With respect to the cluster 
both galaxies have a velocity of ~1000 kms“^ mainly 
towards the cluster center (M87). Thus the atomic gas 
outside the optical radius of NGC 4639 should be af¬ 
fected by ram pressure in the same way as the gas 
of NGC 4654, i.e. it should be mainly located to the 
south-east of the galaxy center. Warmels (1988) ob¬ 
served indeed extended gas at the east of NGC 4639. 
However, this is not observed by Cayatte et al. (1990). 
They found an external emission region in the south¬ 
west of the galaxy center. Moreover Phookun & Mundy 
(1995) did not detect any emission in the east outside 
the optical radius of NGC 4639. Since NGC 4639 is 
located at the edge of the primary beam of Cayatte 
et al. (1990) and Phookun & Mundy (1995) their re¬ 
sults have to be confirmed by observations especially 
dedicated to NGC 4639. 


8. Summary and conclusions 

Numerical models including ram pressure and gravita¬ 
tional interaction are compared to high-sensitivity Hi ob¬ 
servations of the spiral galaxy NGC 4654 in the Virgo 
cluster. This galaxy shows a very extended, low surface 
density tail that does not show clear signs of rotation. 

Four different models are presented, in order to inves¬ 
tigate the origin of this tail: 

— fixed gravitational potential of the galaxy, ram pres¬ 
sure is included, 

— stellar content and dark halo simulated by a non- 
collisional component, ram pressure is included, 

— gravitational interaction, 

— gravitational interaction, ram pressure is included. 
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8.1. Ram pressure as the only perturbation 

The temporal ram pressure profile corresponds to a realis¬ 
tic galaxy orbit within the cluster. Among the simulations 
of Vollmer et al. (2001) with different inclination angles 
i between the disk and the orbital plane, one simulation 
was chosen, which reproduces the observed Hi deficiency 
and the Hi distribution. Limits for these parameters are: 
5000 cm“^(kms“^)^ < Pram < 10000 cm“^(kms“^)^ and 
0° < z < 10°. The timestep of the snapshot compared to 
the observations is roughly given by the projected distance 
of NGC 4654 to the cluster center divided by 1700 kms“^, 
which represents a typical velocity (with respect to the 
cluster mean velocity) of a mildly stripped galaxy in the 
core of the Virgo cluster (Vollmer et al. 2001). With the 
given PA and inclination angle, the azimuthal angle has to 
be chosen such that the radial velocity of the model galaxy 
equals approximately that of NGC 4654. These are very 
severe constraints for the model, reducing the number of 
possible snapshots to a few. 

8.2. Gravitational interaction 

The parameters for the gravitational interaction were cho¬ 
sen such that (i) the observed position and radial velocity 
of NGC 4639 are reproduced and (ii) the observed asym¬ 
metry of the stellar content of NGC 4654 is reproduced. 
Asymmetric spiral arms without disrupting the galaxy can 
only be obtained with a very close, retrograde encounter. 

8.3. Results 

The following results were obtained from the different 
models: 

1. The model that uses ram pressure as the only pertur¬ 
bation shows a tail up to 800 Myr after a strong ram 
pressure event, i.e. after the galaxy’s passage through 
the cluster center. 

2. During the ICM-ISM interaction the non-collisional 
component of the galaxy is heated resulting in a tail of 
low surface density. In the case of a fixed gravitational 
potential the tail has a much higher surface density, 
because of the lower scale height of the disk’s gravita¬ 
tional potential. 

3. An edge-on ICM-ISM interaction can produce an 
asymmetric stellar distribution and thus an asymmet¬ 
ric rotation curve. 

4. An edge-on ICM-ISM interaction cannot produce bar. 

5. The model using a gravitational interaction as the only 
perturbation can reproduce the observed asymmetry of 
the stellar content of NGC 4654. 

6 . The model using a gravitational interaction as the only 
perturbation cannot reproduce the observed, extended 
gas tail of NGC 4654. 

7. Only the model of a mixed perturbation (gravitational 
and ram pressure) can account for all observed proper¬ 
ties (gas and stars) of NGC 4654. In the case of a past 
tidal interaction only a small amount of ram pressure 


is needed to form the observed asymmetries of the gas 
distribution and velocity field. 

8 . Only the comparison with the gas distribution and ve¬ 
locity held can discriminate between the models. 

I thus conclude that NGC 4654 has suffered most prob¬ 
ably a tidal interaction with its companion NGC 4639 
~500 Myr ago. It is now entering the cluster with a ve¬ 
locity of ~I000 kms“^ mainly to the west. It is now 
experiencing a small ram pressure of the order Pram ~ 
200 cm“^kms“^ that is responsible for the observed ex¬ 
tended, low surface density gas tail. 
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Fig. 17. Position-velocity cut through the Hi cube along 
the major axis. Top panel: model of a gravitational in¬ 
teraction alone. Middle panel: model of a gravitational 
interaction together with constant ram pressure. Bottom: 
Hi (Phookun & Mundy 1995). 


Fig. 16. Top panel: model Hi velocity field of a gravita¬ 
tional interaction alone. The contours, in kms“^, are from 
900 (north-west) to 1160 (south-east) in steps of 20 for 
both velocity fields. Middle panel: gas velocity field of the 
model of a gravitational interaction together with constant 
ram pressure with the same contour levels. Bottom panel 
: Hi velocity field of NGC 4654 (Phookun & Mundy 1995) 
with the same contour levels. 
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Fig. 18. Position-velocity cut through the Hi cube along 
the extended tail. Top panel: model of a gravitational in¬ 
teraction alone. Middle panel: model of a gravitational 
interaction together with constant ram pressure. Bottom: 
Hi (Phookun & Mundy 1995). 


Fig. 19. Position-velocity cut through the Hi cube along 
the minor axis. Top panel: model of a gravitational in¬ 
teraction alone. Middle panel: model of a gravitational 
interaction together with constant ram pressure. Bottom: 
Hi (Phookun & Mundy 1995). 





